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•  Army  needs  high-fidelity  multibody  dynamic  modeling  of 
wheeled  and  tracked  vehicles  for  predicting: 

-  Stability  (going  over  a  rough  terrain  at  high-speed,  lane  change,  etc.) 

-  Mobility  analysis  (going  over  bumps  and  potholes). 

-  Fatigue  life  and  durability  of  various  vehicle  components. 
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•  Two  modeling  and  simulation  approaches: 

-  Traditional  solution:  Rigid/Flexible  Multibody  Dynamics  (MBD)  codes. 

-  Alternate  solution:  Explicit  Finite  Element  (FE)  codes. 


•  Advantages  of  Explicit  FE  codes  over  MBD  codes. 

-  Multi-physics  modeling  such  as  fluid-structure  interaction  or  thermal 
effects. 


-  Same  software  environment  to  model  flexible  bodies  and  to  predict 
stresses 

•  However,  accurate  stress  analysis  requires  refined  mesh  and  takes  significantly 
more  CPU  time. 


-  Easily  accommodate  nonlinear  material  characteristics  such  as 
plasticity  and  fracture. 

-  Can  be  easily  parallelized. 


-  Simulation  time  increases  linearly  with  the  DOFs. 
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•  Disadvantages  of  explicit  FE  codes  over  MBD  codes. 

-  Most  explicit  FE  codes  use  an  incremental  updated-Lagrangian 
solution  formulation  which  can  lead  to  solution  drift. 

-  Use  inexpensive  finite  elements  that  use  spurious  modes 
control. 

-  Calculated  joint  reaction  forces  exhibit  high-frequency 
oscillations.  MBD  codes,  on  the  other  hand,  produce  accurate 
joint  forces. 

-  Integration  time  step  (At)  must  be  less  than  a  critical  time  step  for 
stability. 

•  This  requires  a  very  small  time  step  for  stiff  systems.  MBD  codes, 
on  the  other  hand,  can  use  implicit  integrators  that  allow  larger  time 
steps. 
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•  Benchmark  the  two  modeling  and  simulation 
approaches. 

•  Multibody  dynamic  code. 

•  Explicit  finite  element  code. 

•  Benchmark  using  two  multibody  dynamic  systems. 

•  A  7-link  planar  mechanism. 

•  A  spatial  robotic  manipulator. 

•  Comparisons  are  made  of: 

•  Body  motion. 

•  Joint  constraint  forces. 

•  Conservation  of  energy. 

•  CPU  time. 

•  Conclusions  are  drawn  regarding  solution  accuracy  and 
efficiency  of  the  two  codes. 
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Explicit  FE  Code  Formulation 


Equations  of  motion: 


MKx‘Ki  =  F/Ki  +  FarKi 


IkA=t;b  +T.'a  -fe  L 


Integrated  using  trapezoidal  explicit  integration  formula: 


x‘Kj=x‘-;t+0.5A tix'cj+x^) 
xKj  =  XK T  +  0.5  Ar  (x‘Kj  +  x'~f ) 


d'Kj=9'K?  i  0.5  A/  (0'Kj  .^,vi 

A^.=0.5A  t(d'Kj+&KjA') 

R'k  =  J£*  R(A0‘Ki) 
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Explicit  FE  Code  Formulation 


Constraint  modeling. 

-  Penalty  formulation  for  modeling  normal  contact/joint  forces. 


Fc  ■  =  (kP  d  +cP  djdi  )d  f 

FI  *  t 

Cl-  — Xp\  -  Xp 2 • 


where  kp:  penalty  stiffness,  cp\  penalty  damping 


Friction  model. 

-  Asperity-spring  friction  model  is  used  to 
model  joint  and  contact  friction. 
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Multibody  Dynamic  Code 
Formulation 


r1 
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•  Equations  of  motion  : 


•  Constraint  equations: 


M 

d® 

dq 


dq 

0 


ao  . 

o(<7,o  =  o  —q~v 

dq 


•  Forms  a  set  of  Differential-Algebraic  Equations  (DAEs). 


•  Solved  using  implicit  integration  methods  such  as  the 
Backward  Differentiation  Formula. 


-  Can  take  much  larger  time  step  than  explicit  methods. 

-  Very  advantageous  for  stiff  systems. 

-  Time  step  is  only  dictated  by  desired  solution  accuracy. 
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7-Link  Planar  Mechanism 

lvlu  1  V 
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•  2D  mechanism  with  7  rigid  links.  c 

•  1  DOF  system.  / 

•  Driven  by  a  motor  torque  applied  at  point  O.  ^ 


-  Toraue  is  removed  at  time  =  0.1  sec  \ 


Benchmark  Problem  1 : 
7-Link  Planar  Mechanism 


•  Time  step  /  CPU  time 
comparison. 


MBD 

FE 

FE 

Tune  step  (s) 
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0.4E-5 
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CPU  tune  (s) 
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Benchmark  Problem  1 : 
7-Link  Planar  Mechanism 
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Link  K1  rotation  angle  comparison. 


Time  (s) 

Beta  (rad) 

%  Difference 
between  Handbook 
& 

Handbook 

[3] 
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FE 
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FE 
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Benchmark  Problem  1 : 
7-Link  Planar  Mechanism 


9 


MODELING  (mo  SIMULATION,  TESTING  and  VALIDATION 


Link  rotation  angles  comparison 


Link  K1  Rotation  Angle 


August  17-19,  2010 


Link  K3  Rotation  Angle 


Time  (s) 


GVSETS 


UNCLASSIFIED:  Dist  A.  Approved  for  public  release 


Benchmark  Problem  1 : 
7-Link  Planar  Mechanism 


•  Link  velocities  comparison. 


Link  K1  Rotational  Velocity 
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Benchmark  Problem  1 : 
7-Link  Planar  Mechanism 
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•  Joint  force  comparison. 
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Joint  B  X-Force  (N) 


Benchmark  Problem  1 : 
7-Link  Planar  Mechanism 
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•  Effect  of  FE  time  step  on  joint  forces. 


Joint  B  X-Force  on  Link  K3 


Joint  B  Y-Force  on  Link  K3 
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Benchmark  Problem  1 : 
7-Link  Planar  Mechanism 
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•  FE  joint  force  with  ramped  applied  torque  and  spring  pre-load. 


FE:  At=  0.4E-5 
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Potential  Energy  (J) 
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IPJ 


Benchmark  Problem  1 : 
7-Link  Planar  Mechanism 
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Energy  comparison. 

Linkage  Mechanism  Potential  Energy 


Linkage  Mechanism  Energy 
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•  3  rigid  bodies  connected  using  cylindrical  and  revolute  joints. 

•  5  DOF  system. 

•  External  loads:  gravity  and  prescribed  joint  actuator  forces/torques. 

•  End  effector  traces  a  straight  line  with  a  trapezoidal  velocity  profile. 

•  Total  solution  time  =  2.0  sec. 


i  Benchmark  Problem  2: 


Spatial  Robotic  Manipulator 
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•  Time  step  /  CPU  time  comparison. 


MBD 

FE 
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;■  Benchmark  Problem  2: 
Spatial  Robotic  Manipulator 
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Benchmark  Problem  2: 
Spatial  Robotic  Manipulator 
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Benchmark  Problem  2: 
Spatial  Robotic  Manipulator 
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•  Manipulator  Joint  Force  Comparison 
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•  Two  multibody  dynamics  benchmark  problems  were 
solved  using 

-  an  explicit  finite-element  code  and 

-  an  implicit  multibody  dynamics  code. 

•  The  two  codes  predict  the  same  system  motion. 

•  Joint  reaction  forces  predicted  by  FE  code  have  high- 
frequency  oscillations  due  to  the  penalty  method  used. 

•  To  eliminate  high-frequency  force  oscillations  when  using 
FE  code,  applied  forces/moments  must  be  continuous 
and  the  simulation  should  start  from  static  equilibrium. 

•  Implicit  MBD  codes  are  computationally  more  efficient 
than  Explicit  FE  codes  for  stiff  MBD  systems. 
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